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Abstract 



We will present a review of the most interesting results on rare B decays from the B 
Factories, based on the data collected by the BaBar Q] and Belle 2 detectors at asymmetric 
e + e~ colliders at the center of mass energy of the T(4S) resonance. 



1 Introduction 



Rare B decays have always been a standard probe for New Physics (NP) searches as they 
offer a prolific ground where to look for deviations from the SM expectation. The very 
low SM rate of these decays often make them unaccessible with the present experimental 
datasets, unless NP effects enhance the rate up to the current experimental sensitivity. In 
this view, if a suppressed decay is observed, clear sign of NP can be claimed. On the other 
hand, if an upper limit (UL) is set, it can constraint NP scenarios. B Factories provide an 
unique environment where to investigate these processes. The high number of B mesons 
pairs produced often allows to approach the needed experimental sensitivity. Moreover, the 
clean environment and the closed kinematic of the initial state enable to obtain a very pure 
sample where to look for these decays. 

In this work we are going to present a review of results from both BaBar and Belle 
collaboration on the B meson decays B -> h^vP, B+ l+vi , B°J/ip(/> and B+ 
K-ir+ir-/K+K-ir-R 
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2 Analyses Overview 



The analyses presented in this work exploit different reconstruction techniques, depending on 
the particles present in the final state. In the decays B°J/ipcf> and B + — ► K~Tr + ir~ /K + K~ir~ 
all particles decaying from the signal B are detectable, and thus a full kinematic reconstruc- 
tion of the event is possible. On the other hand, B + — > l + vi and B h>*>vv decays possess 
one or more neutrinos in the final state, which clearly can not be detected. This particular 
characteristic calls for different analysis techniques, which allows to deal with the lack of 
information regarding these particles. Typically, the closed kinematic of an e + e~ collision is 
exploited to constrain through energy and four-vector conservation the BB pairs, after both 
particles have been reconstructed. Different approaches can be employed in the selection of 
the B meson which is not decaying into the channel of interest {Btag ) : a totally inclusive re- 
construction is applied on the B tag , without trying to identify its decay products, whenever 
the additional kinematic constraint coming from the two-body nature of the signal B (B S i g ) 
can be exploited, as in B + — » p + vi and B + — » e + vi analyses. The high efficiency obtainable 
with this method has as drawback a poor energy resolution. On the other hand, when more 
than one neutrino is present in the event, a recoil technique is needed: first, the B tag is recon- 
structed in either a semileptonic B s i — > D^*>lu or hadronic Bhad — * DY (Y — n, K) system. 
Then, the channel of interest is searched in the rest of the event (ROE), defined as the set 
of tracks and calorimeter deposits not associated with the B tag . The recoil method allows 
a very high resolution and purity, but has a low efficiency. 

3 Search for b — » svv processes 

3.1 Theoretical Introduction 

In the SM b — > svv transitions occurs through Flavour Changing Neutral Current (FCNC) 
and are therefore forbidden at tree level. As these processes occurs via one-loop box or elec- 
troweak penguin diagrams, they are expected to be highly suppressed. The expected branch- 
ing ratios (B ) are B (B -> K*uv) = (6.8±i;?) x 1CT 6 and B {B -> Kvv) = (4.5±0.7) x 1CT 6 

3. However, this values can be enhanced in NP scenarios, where several mechanisms can 

contribute to the rate. For example, in Ref.il, non _ standard Z o coupling can give rise to 

an enhancement up to a factor 10. Moreover, new sources of missing energy, such as light 

dark matter^ or unparticles if accompanied by a K* , would contribute to the total rate. 

The kinematic of the decay can be described in terms of s„ p = m^ p /m^, where m^p is the 

invariant mass of the neutrinos pai r and tub is the B meson mass. As NP can strongly affect 
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the decay in terms of s„ p shap e 1 1 it is important to not rely on any theoretical model when 
performing the analysis. 

3.2 Analysis of B —> h^vv 

The Belle collaboration analysis exploit the hadronic recoil technique and looks for final states 
with ftM = K + ,K%,K* ,K* + ,Tr + ,n ,p + ,p°,(j), using a 492 ffcr 1 data sample^. B tag candi- 
dates are reconstructed into a charged or neutral £)'*' accompanied by a charged 7r or p or a± 
or mesons are reconstructed as D — > K^tt ,Kgir 7r ,iT!y7r 7r + 7r ,K + -k tt 

and K+tt-tt-tt , while L>W° as £>W° -> K+ir-, K+ir-n , K+ ^ tt_+ tt~ , K%tt° , K%it+^- ,K+ K~ 
and K° s tt-tt + 'k () . The D*-{D ( -*'> ) mesons are reconstructed as D n-(D TT and D°j) and 
the D*+ as D*+ -> D+-/ and D+ -> K%K + and K+R-1T+. No additional charged tracks 
or 7T° candidates are allowed in the event and thus B sig candidates are selected using the 
variable -Eecl = -E-tot — E tec , where E to t and E rec are the total visible energy measured by 
the calorimeter and the measured energy of reconstructed objects including the B tag and the 
signal side candidate, respectively. 

The dominant background comes from BB decays involving a b — > c transition. In order 
to suppress such decays, the momentum of the ft,'*) candidate in the B S i g rest frame P* 
is required to be 1.6 GeV/c < P* < 2.5 GeV/c. The cosine of the missing momentum in 
the lab frame is required to lie between -0.86 and 0.95 to rejects background events where 
particles are missing along the beam pipe. Given the fact that reconstruction efficiencies 



Table 1: Observed number of events N a b s , expected background, ATj,, reconstruction efficiencies e and the UL at 

90% of confidence level for B —* b>*^vv 



Mode 


N obs 


N b 


e(xl0~ 5 ) 




UL 






K* u vi> 


7 


4.2 ± 1.4 


5.1 ± 0.3 


< 


3.4 


X 


10 


-4 


K*+uu 


4 


5.6 ± 1.8 


5.8 ± 0.7 


< 


1.4 


X 


10 


-4 


K+uv 


10 


20.0 ± 4.0 


26.7 ± 2.9 


< 


1.4 


X 


10 


-5 


K°uv 


2 


2.0 ± 0.9 


5.0 ± 0.3 


< 


1.6 


X 


10 


-4 


n + vv 


33 


25.9 ± 3.9 


24.2 ± 2.6 


< 


1.7 


X 


10 


-4 




11 


3.8 ± 1.3 


12.8 ± 0.8 


< 


2.2 


X 


10' 


-4 




21 


11.5 ± 2.3 


8.4 ± 0.5 


< 


4.4 


X 


10 


-4 


p + vv 


15 


17.8 ± 3.2 


8.5 ± 1.1 


< 


1.5 


X 


10 


-1 


<bvv 


1 


1.9 ± 0.9 


9.6 ± 1.4 


< 


5.8 


X 


10' 


-5 



for the UL evaluation are estimated with MC simulation based on the SM, these two cuts 
introduce a SM dependence in the analysis results. The background £"ecl distributions 
are normalized by the number of events in the sideband region. None of the signal modes 
show a significant number of signal events. The observed number of events N b s , expected 
background Nb, reconstruction efficiencies e and th e UL at 90% of confidence level obtained 
with an extension of the Feldman-Cousins method EEO] are shown in Tab. [T] 



3.3 Analysis of B — > K*vv 

The B — ► K*vv from BaBar collaboration is performed in the recoil of both hadronic (HAD) 
and semileptonic (SL) system on a data sample of 413 fb _1 EH. The two different tagging 
strategies provide non overlapping samples whose results can be combined as independent 
measurements. The event selection starts from the B tag reconstruction: in the SL analy- 
sis, neutral D mesons are reconstructed in the K~it + , K~it + it , K~ir + ir~ir + and Kgir + ir~ 
modes. Charged D mesons are reconstructed in the K~ir + ir + and Kgir + final states. In the 
HAD analysis, the -Bhad is reconstructed in i?had — > DY where Y = nir + mK + rKg + qn 
with n + m + r + q < 6 and D is a generic charmed meson. About 1000 different decay chains 
are considered. Charmed mesons are reconstructed in the same final states used in the SL 
analysis, along with the additional channels D + — > K + ir~7r + ir° , Kgir + ir~ir + , Kgir + ir Q . For 
each reconstructed tagging B, a K* is searched in the ROE and reconstructed in the K + tt~ , 
KgTT + or K + tt° mode. Considering that signal events have no additional neutral particles 
produced in association with the K* , one of the most discriminating variable between sig- 
nal and background is, as in the Belle analysis, the extra neutral energy E cxt r&, the sum of 
the energies of the calorimeter neutral clusters not used to reconstruct either the B s i g or 
the B tag ■ In the SL analysis, the signal yield is extracted through a Maximum Likelihood 
(ML) fit to the final -E e xtra distribution, after selection criteria are applied to suppress the 
continuum background. In the HAD analysis, a loose selection is applied and all discrim- 
inants variables (including -E e xtra) are used as inputs for a Neural Network (NN), whose 
output variable NN out is fitted to extract the number of signal events. No significant signal 
is observed in the two analysis and a Bayesian approach is employed to set the UL at 90% 
of confidence level. The UL from the SL analysis are B(B° — » K* Q vv) < 18 xlO -5 and 
B(B+ -> K* + vv) < 9 xl0~ 5 , the UL from the HAD analysis are B(B° -> K* Q vv) < 11 
xlO" 5 and B{B+ -> K*+vv) < 21 xlO" 5 . The combined UL are B(B° -> K*°vD) < 12 
xl0~ 5 , B{B+ -> K*+vv) < 8 xlO" 5 and B{B -> K*vv) < 8 xl0~ 5 . 

These results are at the moment the most restrictive UL on these decay channels and are 
the first completely model independent measurement. 



4 Search for B + — ► l + v processes 

4-1 Theoretical Introduction 

In the SM the purely leptonic B decays B + — > l + v (I = e, fx, r) proceed through the anni- 
hilation of the two quarks in the meson to form a virtual W boson. The branching ratio 
can be cleanly calculated in the SM and is sensitive to the Cabibbo Kobayashi Maskawa 
matrix element V u b and the B decay constant fg, which describes the overlap of the quark 
wave functions within the meson and which is currently the major source of uncertainty in 
the B calculation. Assuming t b = 1.638 ± 0.011 ps, V ab =(4.39 ± 0.33) xl0~ 3 determined 
from inclusive charmless semileptonic B decays^! and fg = 216±22 MeV from lattice QCD 
calculation^, the SM estimate of B(B+ -> r + n ) is (1.59 ± 0.40) x 10~ 4 . Due to helicity 
suppression, B + — > fi + vi and B + — > e + vi are suppressed by factors e /m£ with respect to 
B + — » t + vi , leading to expected branching fractions of B(B + — ► /j, + Vi ) = (5.6 ±0.4) x 10~ 7 
and B(B + — > e + vi ) = (1.3 ± 0.4) x 10 — n . Purely leptonic B decays are sensitive to physics 
beyond the SM, where additional heavy virtual particles contribute to the annihilation pro- 
cesses. Charged Higgs boson effect s m ay greatly enhance or suppress the branching fraction 
in some two-Higgs-doublet models^^l. Moreover, in a SUSY scenario at large tan/3, non- 
standard effects in helicity-suppressed charged current interactions are potentially observ- 
able, being strongly tan /?-dependent and leading to^^l g^+^+^j"^ « (1 — tan 2 /3-^f-) 2 . 
These decays are also potential probes for Lepton Flavour Violation (LFV) in the ratios 
R B /T = B(B+ -» v+vt )/B(B+ -» r+Pt ) and R e J T = B(B+ -» n+vi )/B(B+ -» t+vi )H5I. 

4-2 Analysis of B + — > l + vi in the SL Recoil 

The B + — > Z + z^ analysis from the BaBar collaboration is performed in the recoil of a semilep- 
tonic system on a data sample of 418 fb -1 Q2]. The B tag is reconstructed in the set of semilep- 
tonic B decay mode B~ — > D°l~uX, where / = e, /i and X can be either nothing or a transi- 
tion particle from a higher mass charm state decay, which is not reconstructed (although tags 
consistent with a neutral B decay are vetoed). The D° is reconstructed in the same modes 
as in 13.31 Since r decays before reaching active detector elements, the B + — > t + vi signal 
is searched for in both leptonic and hadronic r decay modes: t + — > e + u e v T , jX^v^Vr , 7r + z? T 
and 7r + 7r P T . The lepton momentum in the B S i g rest frame is used to separate electrons and 
muons from B + — > or _B + — > e + i^ and from r decays. Backgrounds consists primar- 

ily of B + B~ events in which the B tag has been correctly reconstructed and the recoil side 
contains one signal candidate track and additional particles which are not reconstructed by 
the tracking detectors or calorimeter. Typically, these events contain K\ candidates and/or 
neutrinos. In addition, some excess events in data, most likely from two-photon and QED 
processes which are not modeled in the MC simulation, are also seen. Multiple variables are 
used to suppress backgrounds and are combined in two likelihood ratios (LHRs), which are 
probability distributions designed to produce maximum separation between signal and back- 
ground. Among them, t wo o f the most powerful are the ratio of the second and the zeroth 
Fox- Wolfram moment and the cosine of the angle between the B meson momenta and 

the D°l candidate in the T(4S) frame. Also in this analysis, the total energy recorded in the 
detector and not assigned to either the B sig or the B tag E 0Ktia is used to select signal decays. 
The number of expected background events is estimated from the E rmextra sidebands. The 
observed number of events N f, s , expected background Nf,, reconstruction efficiencies e and 
the branching ratio B or UL at 90% of confidence level obtained with the Feldman-Cousins 
method ^ are shown in Tab. [2] 

Combined with the previous measurement in the HAD recoil the B + — > t + vi average 
from BaBar is B(B+ -> t+u) = (1.8 ± 0.6) x 10" 4 . 

4-3 Inclusive Analyses of B + — > l + vi (I = e,/i) 

B + — * l + vi are two-body decays, so the lepton is produced mono-energetic in the B S i g rest 
frame. Thus, in inclusive analyses, the highest momentum lepton in the event is searched 
and assign to the signal side, and all other charged tracks and neutral clusters in the event 



Table 2: Observed number of events N i, s , expected background Nb, reconstruction efficiencies e and B or UL at 

90% of confidence level for B + — > l + vi 



Mode 


N obs 


N b 


e(xl0" 4 ) 




B 






610 


521 ± 31 


10.54 ± 0.41 


(1.8 ± 0.8 


± 0.1) 


ixlO" 4 




11 


15 ± 10 


27.1 ± 1.2 


< 11 x 10" 


" 6 @ 9 


0% CL 




17 


24 ± 11 


36.9 ± 1.5 


< 7.7 x 10 


-6 @ g 


'0% CL 



are used to reconstruct the B tag , without trying to identify the direct decay products of 
the B tag ■ The momentum direction of the reconstructed B tag is used to boost the lepton 
candidate in the B S i g rest frame and thus refine the estimate of the lepton momentum in 
this frame (p*). The two most significant backgrounds are B semileptonic decays involving 
B — > X UtC lv transitions where the endpoint of the lepton spectrum approach that of the 
signal, and non- reso nant qq events. 

Belle analysis^^is on a data sample of 253 fb _1 and applies tight cuts on lepton momen- 
tum in both CM frame (pcm) and B S i g rest frame p* to rem ove B — » X UiC lu backgrounds, 
and exploit the combination of modified Fox- Wolfram moments ^ ^ in a Fisher discriminant 
to suppress continuum events. A cut on the B tag AE = Eb — -Ebcam, where Eb is the recon- 
structed energy of the B tag and -Bbeam the beam energy in the CM frame, is applied to refine 
the selection. The final selection efficiency is (2.2 ± 0.1)% for B + — > and (2.4 ± 0.1)% 
for B + — » e + vi . The yields are extracted from a ML fit to the beam-energy constrained 
mass M oc = y/ E£ — p% distributions, where ps is the reconstructed momentum of the 
B tag , and 4.1 ± 3.1 signal events are observed for the muon mode and -1.8 ± 3.3 for the 
electron mode. The 90% confidence level for the upper limit on the branching fraction Bgo 
are defined by 0.9 = f® 90 C{B)dB/ J °° C{B)dB and are B (B+ -> ^v x )< 1.7 x lO" 6 and B 
(B+ -> e+vi ) < 9.8 x 10~ 7 

The 90% confidence level UL for the electron mode is currently the most stringent mea- 
surement available . 

BaBar analysis^] uses an integrated luminosity of 426 fb _1 and choose to combine five 
different topological and kinematical variables, optimized separately for each mode, in a 
Fisher discriminant for qq background suppression. A cut on the B tag AE is applied for 
the muon mode in order to remove continuum background, while two linear combination of 
Btag AE and B tag Pt are used for the electron mode in order to reject also the background 
coming from two-photon events. The final selection efficiencies are (6.1 ± 0.2)% for the 
muon mode and (4.7 ± 0.3)% for the electron mode. The two-body nature of the decay is 
exploited by combining p* and pcu in a second Fisher discriminant, whose output pfit is 
used in combination with the B tag rriES — V-^bcam — \pb \ 2 a ML fit to extract the yields. 
The number of observed signal events is 1.4 ± 17.2 for B + — > [i + vi and 17.9 ± 17.6 for 
B + — » e + vi . The 90% confidence level ULs are evaluated with a Bayesian approach and are 
B (B+ -> fi+vi )< 1.0 x 10~ 6 and B (B+ -> e+vi )< 1.9 x 10" 6 

The 90% confidence level UL for the muon mode is more restrictive than any previous 
measurements. 

5 Search for B° — > J/ip(f> decay 
5.1 Theoretical Introduction 

Studies of exclu sive B meson decays to charmonium play an important role in explor- 
ing CP violationE2l23] anc J establish the Cabibbo-Kobayashi-Maskawa picture. The decay 
B° — > J/ip4> is expected to proceed mainly via a Cabibbo-suppressed and a color-suppressed 
transition (b — > ccd) with rescattering. In B decays, effects p resuma bly due to rescattering 
have been seen in various decay processes, as B° — > D~K + l ^4 | 25 |^ anc j th e y ma y pi a y an 
important rol e in understanding patterns of CP asymmetries in B decays to two charmless 
pseudoscalars^Sl. 



5.2 Analysis of B° — » J/tpcj) 

The Belle analysis of B° — > J/ip(f> decay is performed on a data sample of 605 fb _1 ^l 
Candidates for B Q — > J/ipcj) decays are reconstructed from the decay J/tp — » = e,/i) 

and <fi — > if + i4T~. £?° are selected through the usual Mf, c and AE variables and the signal 
yield is extracted through a ML fit to the AE distribution. 

The dominant background comes from BB events with a B decays to a J/ip, in particular 
B° -> J/ipK*° (892) [-> and -> J/r/>.Ki(1270)[-> if-Tr+Tr / - ]. In both cases, a 

pion is misidentified as a kaon, and in the latter case the other pion is missed. The former 
has a peak at AE ~ 0.1 GeV, while the latter has a broad peak in the negative AE region. 
These background are taken into account into the ML fit, as well as the remaining not- 
peaking combinatorial background. The number of observed signal events is 4.6lJ'5 with a 
significance of 2.3 a and the 90% confidence UL, extracted by a frequentistic method using 
ensembles of pseudo-experiments, is B(B° — > J/ip<fi) < 9.4 xl0~ 7 . The UL is the most 
restrictive up to date and improve of about one order of magnitude the previous result^. 

6 Search for b — > qqd/qqs Processes 

6.1 Theoretical introduction 

b — ► qqd/s transitions are highly suppressed in the SM: compared with the penguin (loop) 
transition b — > qqd/s, they are additionally suppressed by the small Cabibbo-Kobayashi- 
Maskawa matrix element factor | VtdV t * | ~ 3 x 10~ 4 , l eading to a predicted branching fractions 
of only O(10~ 14 ) and O(10~ n ), respectively ^ ! *^ . These branching ratios can be signifi- 
cantly enhanced in SM extensions as the minimal supersymmetric standard model (MSSM) 
with or without conserved R parity, or models containing extra U(l) gauge bosons'^. Ob- 
servation of the decays B~ — » K + ir~ir~ and B~ — > K~K~ir + would be clear experimental 
signature for b — > dds and b — > ssd quark transitions, which have already been searched in 
these and other decay modes without success yet. 

6.2 Analysis of B~ — > K + tt~tt~ / K~ K~ir + 

The BaBar B~ — > K + tt~tt~ / K~ K ~ir + analysis uses 426 fb _1 of integrated luminosity^]. 
B~ —> K + ir~ir~ (K~ K~tt + ) candidates are selected by combining a charged kaon (pion) 
candidate with two charged pion (kaon), each of which has charge opposite to the kaon 
(pion). In order to avoid potentially large source of background arising from B decays 
mediated by the favored b — ► c transitions, B decays in which the pairs of daughter tracks 
have invariant mass combinations in the range 1.76 < < 1.94 GeV/c 2 , 2.85 < m K7r < 

3.25 GeV/c 2 and 3.65 < < 3.75 GeV/c 2 . These vetoes remove events containing the 

decays D° — > K~ir + , J/tp — + l + l~ and <fi(2S) — > l + l~, respectively, where the leptons in the 
J/ip and <p(2S) decays are misidentified as pions and kaons. The final selection efficiencies 
are 21.6% for B~ — > K~ir + Tr + and 17.8% for B~ — > K ~ K~-k + . Continuum events represent 
the dominant background: in order to discriminate against it, five variables, comprising 
topological quantities as well as the flavour properties of the recoiling B meson and the 
proper time difference between the two Bs, are combined into a NN. Residual backgrounds, 
arising from decays topologically similar to the signal but with some misreconstruction, can 
be divided in five categories, depending on their shape in tues an d AE, and are taken into 
account in the fit to the yield. The number of signal events is extracted through a ML fit 
to m ES , AE and the NN output NN out and is found to be 22 ± 43 for B~ -> K^tt + t: + 
and -26 ± 19 for B~ — > K~ K~n + . A frequentistic Feldman-Cousins method is employed 
to obtain the 90% confidence level UL, which are B(B~ — > K~tt + tt + ) < 7.4 xlO -7 and 
B(B~ — > K~ K ~7r + ) < 4.2 x 10~ 7 and improve of about a factor 3 previous measurement^^. 

7 Conclusions 

In this work we have presented a review of the most interesting results on rare B decays at 
the B Factories from BaBar and Belle experiments. 



Rare decays posses high interest, as they are standard probe for NP searches, given the 
low decay rate expected in the SM. Their study is complementary to the direct exploration 
of the energy frontier and in some cases can access even higher energy scales. We have 
seen how the improved analysis techniques and the huge integrated luminosity from both 
BaBar and Belle experiments allow today to reach 0(1CP 6 — 1CP 7 ) sensitivity and how, even 
if only UL, the results on these decay are already able to pose interesting constraints on 
various NP scenarios. Nonetheless, decays with undetectable particles in the final state will 
not be measurable at the LHC and a Super Flavour Factory will be needed in order to obtain 
improved measurements. 
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